Macroscopic superpositions of superfluid flows 
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We present a scheme for creating macroscopic superpositions of the direction of superfluid flow 
around a loop. Using the Bose-Hubbard model we study an array of Bose-Einstein condensates 
trapped in optical potentials and coupled to one another to form a ring. By rotating the ring so 
that each particle acquires on average half a quantum of superfluid flow, it is possible to create a 
multiparticle superposition of all the particles rotating and all the particles stationary. Under certain 
conditions it is possible to scale up the number of particles to form a macroscopic superposition. The 
simplicity of the model has allowed us to study macroscopic superpositions at an atomic level for 
different variables. Here we concentrate on the tunnelling strength between the potentials. Further 
investigation remains important, because it could lead us to making an ultra-precise quantum-limited 
gyroscope. 
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Superpositions are one of the defining differences be- 
tween classical and quantum mechanics. To test whether 
quantum mechanics can describe the macroscopic world, 
which would normally be described classically, we will 
look for superpositions in larger systems. Multiparti- 
cle superposition states have been observed in a num- 
ber of systems including photons [H, Ceo molecules 0, 
and the internal state of four 9 Be + ions Experi- 
mental signatures of larger scale quantum phenomenon 
were shown when Rouse et al. Q observed resonant tun- 
nelling between two macroscopically distinct states in a 
superconducting quantum interference device (SQUID). 
The observed tunnelling was between states of different 
flux or opposite currents flowing around a loop. Macro- 
scopic systems consist of approximately 10 10 particles or 
have a macroscopic measurable quantity associated with 
them. The currents measured in the SQUID consisted 
of approximately 10 9 Cooper pairs and produce a mea- 
surable magnetic flux, meaning tunnelling between two 
macroscopically distinct states had been achieved. Sim- 
ilar systems have also been used to show cat states can 
be made 0] . 

Bose-Einstein condensates (BECs) are a promising sys- 
tem for realising similar results. They are composed of 
10 3 — 10 7 atoms with a high proportion in the same quan- 
tum state and are sufficiently cold to undergo a quantum 
phase transition from superfluid to Mott insulator [7|. 
They also have significant advantages over SQUIDs since 
they are highly controllable: the coupling between con- 
densates and the strength of the interactions between 
atoms can be tuned over many orders of magnitude, there 
are few imperfections and near perfect lattices can be cre- 
ated. This enables us to develop a simple model to inves- 
tigate macroscopic quantum effects SQUIDS have the 
advantage that a precise magnetic field can be applied to 
the system, which produces an easily controllable phase 
around the loop. This maybe more difficult for BECs. 

There have already been a number of theoretical pro- 



posals for producing cat states with BECs in a range of 
different set-ups Q. In this paper, we present a scheme 
for producing a multiparticle superposition of different 
superfluid flow states in a ring of coupled BECs. It en- 
ables us to study quantum effects on an atomic level us- 
ing the Bose-Hubbard model Q rather than the macro- 
scopic approach used to describe the SQUID H @] . We 
show that the macroscopic superposition can be made 
for larger numbers of atoms by changing the tunnelling 
between the sites. Cat states are important because it 
provides a direct manifestation of quantum mechanics at 
the macroscopic level in a new system. As discussed by 
Leggett such states are important for testing the 
limits of validity of quantum mechanics. Once the sys- 
tem is understood it maybe possible to design a device 
that can make quantum-limited measurements of angular 
momentum or, equivalently, ultra-precise gyroscopes. 

Our system consists of condensed atoms trapped in an 
optical potential of three sites in a ring configuration. 
Each site is coupled to its neighbours by quantum me- 
chanical tunnelling through the potential barriers sepa- 
rating them. By Larmor's theorem ll|, the analogue of 
applying an external flux to create a superposition in the 
SQUID experiments is to rotate the ring, which is equiv- 
alent to applying a phase around the ring. The same 
effect could be achieved by rotating the optical poten- 
tial directly, producing a flow of atoms round the loop 
by using Bragg scattering to imprint phases on the lat- 
tice sites [12l . 1 1 31 ] , or producing an effective magnetic field 
using two resonant laser beams 13, Hj]. 

We can write the spatially dependent wave function for 
the condensed atoms as ipo (x) = e 1 **^ | ipo (x) | , where $ (x) 
is the phase of the condensate at position x. When the 
phase is not constant throughout the condensate there is 
a velocity field associated with it, v(x) = (h/m)V<fr(x). 
The phase of the condensate must be uniquely defined at 
all points round the loop. For a linearly varying phase, 
corresponding to a flow state, this restricts the phase 
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variation around the loop to integer multiples of 2ir. In 
this paper we shall be looking at situations beyond the 
mean field theory where the atoms are in a superposition 
of flow states. These are found by full diagonalisation of 
the Bose-Hubbard Hamiltonian. 

The applied phase between the sites can be incorpo- 
rated into the Bose-Hubbard Hamiltonian Q by includ- 
ing appropriate phase factors in the coupling terms. This 
gives the 'twisted' Hamiltonian 

H = - J[e^ (a f 6 + tfc + c+ a) + (tfa + Jb + a^c)} 
+U(a^a 2 + b^b 2 + c^c 2 ), (1) 

where a, b and c are the annihilation operators of atoms 
in the three sites, U is the on-site interatomic interac- 
tion strength and J is the tunnelling strength between 
adjacent sites. The phase factors e ±l * in the coupling- 
terms are known as Peierl's phase factors. We note that 
cj) does not have to obey the phase matching condition 
because it represents the sites and not the condensate, 
and it is related to the angular momentum of the sites, 
L, by </> = 2nL/3h. 

It is convenient to consider a new orthogonal basis of 
operators {a,(3,j} that incorporates the phase match- 
ing condition. One way to write this is in the quasi- 
momentum basis, or flow basis, 

a = (a + b + c)/y/S, 

f3 = (a + be 1277 ' 3 +ce l4 "/ 3 )/V3, 

7 = (a + be- 12 *' 3 + ce- l4 " /3 )/V3. (2) 

The new basis respectively correspond to annihilation of 
an atom with zero flow, one quantum of clockwise flow 
and one quantum of anticlockwise flow. They follow the 
usual commutation relations and the quasi-momentum 
conservation rules for a periodic potential. Throughout 
this paper we use the convention that a positive phase 
variation corresponds to clockwise flow. Using this we 
can rewrite the 'twisted' Hamiltonian in flow representa- 
tion, 

H = —J{ (2a t a -ft/3- 7 f 7 ) cos <p + V3(/3 f /3 - 7*7) sin 0} 

+ ^{Q,t2 a 2 + pXIpl + 7 f2 7 2 + 4(^0^ + a t a7 t 7 

+/3W7) + 2(a 2 /3 t 7 t + /3 2 aV + 7 2 a t /? t + h.c.}. (3) 

From Eq. J^J we see that the eigenstates of the system, 
when J7/J<C 1, are just flow states. Even in this limit 
Eq. [3] is still valid, because the atoms still tunnel through 
the barriers [8[. For f7/J«l, there is coupling between 
the different flow states and when U/J^>1 the system is 
in the Mott regime where each site acquires exactly the 
same number of atoms. 

Our procedure for creating cat states is as follows. 
Condensed atoms are trapped in the optical ring de- 
scribed above and the tunnelling strength between the 
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FIG. 1: The 4 lowest energy levels of the system for three 
atoms are shown as a function of phase, <f>, between two sites, 
in a ring of three, is varied, and U/J = 0.5. A clear anti- 
crossing between the lowest two energy levels is seen. 

sites is adiabatically reduced by increasing the intensity 
of the trapping light. This results in squeezing of the 
number of atoms at each site until eventually a Mott 
transition takes place 0]. This makes the system very 
stable to flow state excitations when rotated as can be 
seen by the absence of any phase dependence in Eq. [3] 
when J = 0. Next, the loop is rotated at a rate cor- 
responding to half a quantum of flow per particle and, 
while it continues to be rotated, the tunnelling through 
the barriers is adiabatically increased. As long as the 
system evolves adiabatically the system will stay in the 
ground state, so the experiment can be considered at zero 
temperature. We claim this procedure creates a cat state 
of the form, 

(|7V,0,0) + |0,7V,0))/V2, (4) 

where the three terms in the ket represent the number of 
atoms in the a, (3 and 7 flow states respectively. 

Once the final state has been produced we need to 
show it is a cat state. One signature of an entangled 
state is an anti-crossing of the two lowest energy levels. 
In the SQUID experiment of Friedman et al. Q , an anti- 
crossing of two levels was observed. This signature can 
be understood in the context of work on entanglement 
witnesses [l7j . The lowest possible energy for separa- 
ble states can be found using a variational approach and 
gives a crossing of the two lowest energy levels when the 
ring is rotated at a speed corresponding to half a quan- 
tum of flow per particle. This is because the zero flow 
state and the rotating flow state both have the same flow 
speed relative to the potential and therefore have the 
same energy. If the ground state has a lower energy than 
the lowest possible for separable states, we are forced to 
conclude that the state is entangled. 

The energy eigenvalues of the resultant states are plot- 
ted against phase between two sites in Fig. [1] for three 
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atoms, with final values of UJ J = 0.5. The probability of 
being in flow states \N a , Np, AT 7 ) was calculated, where 
N a , Np and iV 7 are the number of atoms in the a, j3 and 
7 flow states respectively. At <f> — the ground state has 
a probability of 0.99 of being in state |3,0,0). The devi- 
ation from 1 is due to interactions between the particles. 
For increasing values of phase there is only a slight drop 
in the probability of this state until <fi 7r/3. At this 
phase the number of atoms in the system becomes im- 
portant. For commensurate numbers of atoms (i.e. the 
ratio of atoms to lattice sites is an integer) there is an 
anti-crossing of the two lowest lying energy levels, which 
is seen in Fig. [T] We have shown that the energy gap 
persists over all values of J when U 0. This means 
that it is always possible to evolve the system sufficiently 
slowly to remain in the ground state. For {7 = the an- 
ticrossing is not observed and the ground state becomes 
degenerate with other flow states at <fr = Tr/3. For exam- 
ple, if the system contained three atoms there would be 
a 4 fold degeneracy of the |3,0,0), |0,3,0), |2,1,0) and 
1 1, 2, 0} flow states. For systems with interactions, as in 
the case shown in Fig. [TJ the energy of states with atoms 
of differing flow become raised. This allows a metastable 
superfluid flow to form in the ground state, which is 
needed to create a cat state In the case of Fig. [TJ the 
probability that the ground state is in the superposition 
(|3, 0,0) + |0,3,0))/a/2 at = 7r/3 is 0.93 and the exited 
state is in the superposition (|3,0, 0) — |0, 3, 0) )/ a/2 is 
close to 1. This shows a good cat state has been formed. 
If the phase is increased further the ground state becomes 
a state close to |0, 3, 0). Non-commensurate numbers pro- 
duce no anti-crossing and there is just a rapid change of 
the ground state from the state close to | TV, 0,0) to the 
state close to |0, N, 0). We shall see shortly how it is pos- 
sible to overcome this difference between the commensu- 
rate and non-commensurate cases. The other two lines 
in Fig. Q] are states close to |2, 1,0) and |1,2,0). Sim- 
ilar results for larger numbers have been obtained, but 
to get similar probabilities of superpositions U must be 
reduced to avoid significant amplitudes in higher energy 
states. For 30 atoms to have the ground state with a 
probability greater than 0.9 of being in a superposition 
U/J must be less than 0.1. 

There is already a substantial amount of work done on 
similar systems at the mean field level. Phenomena such 
as 'swallow tails' can be observed in the energy levels in 
figures similar to Fig. Q] H EH ■ This is not observed 
in our system for any chosen parameters and suggests 
a possible limitation of the mean-field approach. Direct 
comparisons with mean field should await a fully corre- 
lated calculation for larger systems. 

The flow distribution for 30 atoms has been plotted 
in Fig. [2] (a) in the Mott regime with = 7r/3. The 
distribution equally favours the states 1 30, 0, 0), |0, 30, 0) 
and |0, 0, 30) and has a preference of atoms being in the 
same state, which is shown by the largest peaks being in 





FIG. 2: (a) shows the ground state of the 'twisted' Hamil- 
tonian for 30 atoms, where U/J = 1000 and (f> — tt/3. (b) 
shows the ground state for 30 atoms, where U/J = 0.1 and 
4> = 7r/3. The two large peaks demonstrate a cat state. Each 
bar represents the probability of finding the ground state in 



flow state \N- 
of atoms. 



- Np — 7V 7 , N/3, iV 7 ), where N is the total number 



the corners. We will show how this state can be used to 
create a cat. For non-commensurate number of atoms the 
probability distribution favours atoms being in |7V, 0, 0) 
(or |0, N, 0) as there is no preference at </> = 7r/3), so only 
peaks at one corner and is zero in the others. This state 
is not a suitable starting point for creating a cat by the 
process described here. 

When the tunnelling is increased, so the ratio U / J < 
0.1, for the commensurate case, the system is allowed to 
go into states \N, 0,0) and |0, TV, 0) , because there was 
an amplitude in both these states, when the tunnelling 
was weak. This gives rise to the superposition shown in 
Fig, [2] (b), where a significant proportion of the ampli- 
tude takes the form given by Eq. [4j However, for the 
non-commensurate case, when {//</< 0.1, there is no ini- 
tial probability of being in |0, TV, 0) (or | AT, 0,0)), so the 
system is forced into |iV, 0, 0} (or |0, AT, 0)) and no cat 
state is formed. A graph would show only one bar at the 
point representing \N, 0,0) (or 1 0, iV, 0} ) of height close 
to 1. 

A simple way of overcoming the difference between the 
results for commensurate and non-commensurate num- 
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bers of atoms is to have slightly different tunnelling be- 
tween the sites. As you can see from Eq. [3] the cou- 
pling between different flow states is due to the on site 
interaction term. All terms in Eq. [3] conserve flow (or 
quasi-momentum) and only |7V, 0, 0) and |0, iV, 0) have a 
coupling path through intermediate flow states for com- 
mensurate numbers. For example, the state |3, 0,0) can 
couple to |0, 3,0) through |1,1, 1) by applying a 2 ftrf 
then aft 7. Without a coupling path there is no way one 
flow state can evolve into another flow state, so no super- 
position is possible between these states. With different 
tunnelling between the sites, coupling between flow states 
where only one atom changes momentum is allowed, so 
\N, 0, 0) and |0, N, 0) can now couple through intermedi- 
ate flow states. For example, the state |2, 0, 0) can couple 
to |0, 2, 0) through |1, f , 0), which was previously not al- 
lowed. The strength of the coupling depends on, among 
other things, the tunnelling strength, which ultimately 
depends on the number of atoms. In an experimental 
realisation such imperfections will inevitably be present. 
Now when the tunnelling strength is increased, the flow 
distributions peak at \N, 0, 0) and |0, N, 0), so a superpo- 
sition similar to that shown in Fig. [2] (b) is observed. The 
ground and excited states are close to the superposition 
states |iV,0,0)+e l >,7V,0) and \N, 0, 0) +e i ( 9 +^\0, N, 0) 
respectively, where the phase factor depends on the num- 
ber of atoms in the system. The uneven tunnelling, there- 
fore, aids greatly the cat making process, because there 
will be little distinction between commensurate and non- 
commensurate numbers of atoms in realistic situations 
and we should create cats for all numbers. A dynamical 
simulation of the increasing of the tunnelling strength, 
while the loop is rotating at a rate corresponding to half 
a quantum of flow, is shown in Fig. [3] for six atoms. 

The effect of number and tunnelling strengths can be 
interpreted in a more intuitive way. For U S> J, and 
equal values of J around the ring, the ground state of 
the system is highly number squeezed. For the commen- 
surate case, there is no net flow due to the high energy 
cost for doing so. This is why Fig. [5] (a) is symmetric. 
In the non-commensurate case, there are one or two ex- 
tra atoms that can freely tunnel through the barriers. 
These particles have a net momentum and so give rise to 
a skewed momentum distribution. For example, if there 
are 31 atoms, the system can be written as, 

Mott Superfluid 

(aW) 10 (a f + e l V + e 2l V) |0, 0, 0, ), (5) 

where the first term creates the 30 atoms in the Mott 
state and the second term creates the one superfluid atom 
that freely tunnels through the barriers. The value of 8 
determines the momentum state of the extra particle, 
which must satisfy phase matching. 

When the tunnelling strengths are unequal the dis- 
tribution for commensurate atoms is hardly changed. 
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FIG. 3: Shows the probability that the ground state is in the 
flow superposition (|iV, 0,0} + |0, iV, 0))/v2 as a function of 
J/U for a total of six atoms. 

However, for the non-commensurate case, the situation 
changes dramatically. The atoms not trapped by the 
potential see a potential barrier and can either tunnel 
through it or be reflected from it. The condition of 
phase matching means that both these outcomes are both 
equally favourable when the phase across each pair of 
sites is <j) f=a 7r/3 and so a superposition results. For ex- 
ample, for 31 atoms, the state is written, 

(aW) 1( V + 6*^)10,0,0), (6) 

where tp is some phase. As the potential barriers are 
lowered more atoms become free and also see the barrier, 
so contribute to the superposition. 

There are several factors that determine how good a 
cat state can be produce and how accurately we need to 
control the rate of rotation of the ring to produce a cat. 
A discussion of how the tunnelling strength and number 
of atoms effects the cat has been presented in this paper. 
Future work will look at how the tunnelling strength and 
the number of atoms effects the required accuracy of the 
phase to produce a good superposition and whether long 
range interactions will produce more stable cats. The 
system only considered three points that atoms could oc- 
cupy, which is very different to a SQUID that has one or 
several regions that are superfluid. It maybe interesting 
to develop a BEC system that is similar to this, because it 
could prove to be more stable than the three site system. 
Measuring whether a cat state has been formed cannot 
be done directly, however, the anti-crossing could be ex- 
perimentally measured to give evidence of a cat using 
spectroscopic techniques. This is analogous to the read- 
out employed in the SQUID experiments and a detailed 
study of this will form the subject of future work. 

In conclusion, we have demonstrated a scheme for cre- 
ating macroscopic superposition states in the direction of 
the superfluid flow of Bosc-Einstcin condensates around 
a loop. This straightforward scheme relies on being able 
to raise, lower, and rotate an optical potential. All these 
techniques are within reach of current experiments. We 
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have seen that if the tunnelling between each pair of sites 
is exactly equal, cats states are created only for commen- 
surate numbers of particles. However, any slight differ- 
ences in these tunnelling strengths, as we would expect 
to arise naturally in an experiment, means that cats can 
be created for any number of particles. Furthermore, the 
appearance of an anti-crossing between the two lowest 
energy levels provides a clear experimentally accessible 
signature that an entangled state has been formed. The 
simplicity of the model has allowed us to study cat states 
at an atomic level and future work will allow us to inves- 
tigate why cat states are hard to make. This system is 
of great interest as it may have significant technological 
applications especially in the field of metrology and give 
new insight into quantum mechanics. 

This work was supported by the United Kingdom EP- 
SRC and by the Royal Society and Wolfson Foundation. 
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